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Abstract 

A new scheme of the iterative perturbation theory is proposed for the strongly correlated elec- 
tron systems with orbital degeneracy. The method is based on the modified self-energy of Yeyati, 
et al. which interpolates between the weak and the strong correlation limits, but a much simpler 
scheme is proposed which is useful in the case of the strong correlation with orbital degeneracy. It 
will be also useful in the study of the electronic structures combined with the band calculations. 



In the study of strongly correlated electron systems, the dynamical mean-field theory 
(DMFT)Q is recognised to be a very powerful tool. In DMFT, the lattice proplem is reduced 
to a single magnetic impurity embedded in an effective medium which should be determined self- 
consistently. This impurity problem can be solved, e.g., by the numerial renormalization group 
method, [gj which however needs a rather heavy numerical calculation. Therefore, the iterative 
perturbation theory (IPT)|j| is used frequently. This method adopts the modified perturbation 
theory (MPT),§] which uses the second-order perturbation for the calculation of the self-energy of 
the correlated electrons but modifies it so as to interpolates between the weak and the strong corre- 
lation limits in the electron-hole symmetric and asymmetric cases phenomenologically. Despite its 
approximate nature, the method can successfully decsribe both the metallic and the Mott-insulator 
states for the single-band Hubbard model. || 

In the heavy fermion materials, [Q it is of much importance to take account of the orbital 
degenracy. For example, the 4f state with J = 5/2 in the cerium compounds has six- fold degeneracy 
if the crystalline field splitting can be neglected. The Fermi liquid theory is constructed for such 
systems by Yamada, et al.,[[| and a numerical calculation via the self-consistent second-order 
perturbation theory (SCSOPT) has been performed by Kontani and YamadaQ for the periodic 
Anderson model (PAM) with the degeneracy N = 6. However, this calculation is applicable only 
to the case of weak correlation. 

Yeyati, et al.|| proposed an extension of MPT to the case of multi-orbitals and applied it to the 
quantum dot problem, which is equivalently described by the impurity Anderson model with the 
orbital degeneracy. Their method reproduces the atomic (strong correlation) limit for any number 
of localised electrons, but is rather complicated even for N — 4 (spin degeneracy times the double 
degeneracy of orbitals). In the application to the heavy Fermion systems, the Coulomb repulsion 
U between f electrons is usually very large. Therefore, it will be useful if one can find out a simpler 
interpolation scheme than Yeyati's for such a case. 



In this Letter, I propose a new interpolation scheme for the self-energy by taking account of 
only f° and f 1 states of f-electrons in the ground state. It interpolates between the weak and the 
strong correlation limits within the above Hilbert space. The transition to f 2 state is included as 
an excitation. The Fermi liquid properties are completely realised at low temperatures and low 
energies. Therefore, it would be very useful if it is combined with the density of states obtained 
from the band calculations. 

Following Kontani and YamadaQ, we start from PAM with the iV-fold degeneracy for f- 
electrons written as 

H = Yl ekC k<T C k<r + Eff+J lm + U ^ n S^ n S^m 
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Here, Nq denotes the number of lattice sites, each of which is labeled by the index i, and Ef, e k 
and Vko-m are the f-level, the conduction electron energy and the mixing, respectively. The indices 
I and m denote the z component of the angular momentum of the f-states and a the spin of the 
conduction electrons. For 4f electrons with the total angular momentum J = 5/2, V karn is given 

by 

v kam = V^v^i-^J IZ ^S liim _i a Y£ i V*,<M (2) 

if the conduction electrons are assumed to be in the plane wave states. V denotes the strength 
of the mixing. In this Letter, the paramagnetic state with respect to spin and orbital is assumed. 
The on-site Green's funciton for f-electrons Gf(e) = ^2 k G /(k, e) can be written in the form 

(N-2)/N 1 „ 2JN 

°' w " «-*,-=,(«) + 35? (t) _^- (3) 
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The k-dependence of the self-energy is neglected by using the local approximation, which is justified 
in the limit of large spatial dimensions d — > oo.j|] When the self-energy is calculated by the second- 
order perturbation theory (SOPT), S/(e) is given by %(e) = (N - l)Un f + E^ 2) (e), where E^ 2) (e) 
is calculated at zero temperature as 

= {N ^T 2 J de'de"G f (e')G f (e")G f (e + e'- e"), (4) 

and n/ denotes the number of f-electrons per site and orbital. In SOPT, the unperturbed (Hartree- 
Fock) Green's function is substituted to G/(e)'s. However, the Coulomb interaction between 
f-electrons is so large that SOPT is meaningless. Therefore, the self-consistent SOPT (SCSOPT) 
is used in ref. |?| where the self-consistently renormalised Green's function was used for G / (e) 's in 
eq.Q. It is well known, however, the density of states in SCSOPT can not reproduce the lower and 
upper Hubbard bands for N — 2 when the correlation is large. |)| In Fig. 1, I show my calculation 
of the f-electron density of states p/(e) = — 7r _1 ImG/(e + i0 + ) for PAM. The shape of the density 
of states of the conduction band is assumed as a semicircle of the center at Eq and the width W, 
which is taken as the unit of energy. The calculation is done for W = 1, V = 0.25, Ef — 0.15, 
Eq = 0.5 and U = 0.3. A small but finite imaginary part is introduced in the denominator of 
the Green's function. The number of f electrons n/ is obtained as 0.101 per orbital. The mixing 
between the f and the conduction electrons opens the mixing gap, within which one can see the 
sharp peak that comes from the first term in eq.(||) representing the unmixed part of the A^-fold 
degenerate f states. The separation of the lower and the upper Hubbard peaks of the density of 
states are not seen here, which appears only when the iterative perturbation theory (IPT) is used 
for larger U (see below). 

In DMFT, the magnetic atom in the central cite can be regarded as embedded in the effective 
medium expressed by the cavity Green's function G/(e), which is related to the on-site Green's 
function G/(e) asjj] 

G / (c) = (G / ( e )- 1 +E / (e))- 1 . (5) 



IPT calculates the self-energy of an f-electron due to the Coulomb interaction on this impurity site 
by SOPT using Gf(e) as an unperturbed Green's function. This works well only for the symmetric 
case, where even the Mott transition is properly described. || In the asymmetric case, however, it 
fails to reproduce the correct atomic limit. Yeyati, et al.|| proposed to modify the second-order 
self-energy in a phenomenological manner, by which it becomes possible to interpolate between 
the weak and the strong correlation limits. When N = 2, the following form is assumed 

The parameters A and B are determined so as to reproduce the correct |e| — > oo and U — > oo 
limits, respectively: 

A =^-^4 (7) 
n/(l-n/) 

_ U(l-n f )+E f -E f 

where rif — n[Gf] and fif — n[Gf\. The functional n[G] means the electron number calculated by 

n[G] = [ de/(e) (-- J ImG(e + iO+). (9) 



Ej = Ef + Urif denotes the Hartree-Fock level. We call this scheme as the modified perturbation 
theory (MPT) . Recently, I discussed its microscopic basis in |l4j . 

Since the modified form of the self-energy is not conserving, it does not satisfy the Luttinger 
sum rule go), 0: 

Here fi (= 0) denotes the chemical potential and 6(x) denotes the step function. _Ek- denotes the 
lower branch of the quasi-particle band and I have assumed the case that /z lies in the lower band. 
n c denotes the number of conduction electrons. It is possible to approximately satisfy the sum rule 
by subtracting the zero-energy part as|ll| 

S>(e) = E / (e)-E / (0). (11) 

We follow this prescription in the present Letter. If the sum rule must be satisfied explicitly, one 
can regard Ef as a free parameter and determine it by the Luttinger sum rule.J^, ^] 
For N > 2, the f-electron Green's function in the atomic limit is given byj^] 

Gm(£) " 7=E} e-E f -U + "' 

, <iL(#m)»/) (12) 



e-Ef- (N-1)U' 

where the suffices to, £, s, • • • label the orbital states of f-electrons. Yeyati, et al.[|| found a scheme 
which interpolates between this G^ n {e) and the second-order perturbation theory. Their formula, 
however, needs the average values of the products of the f-electron number operators in the different 
orbitals, (n m n£n s • • •), which must be determined self-consistently. The procedure is complicated 
enough even for N = 4. 

In the realistic heavy Fermion materials, it is quite often the case that U is large enough 
compared to the resonance width of the f-states. Then the ground state is restricted to f° and f 1 
states in the case of Ce, so that the above Green's function can be approximated by 

G ^ e >~ J=E t + e-E f -U' {3) 



since the terms which include more than two rig's vanish. The first and the second terms in eq. 
represent the transitions f 1 — > f° and f 1 — > f 2 , respectively. 

In the same manner as TV = 2, we start from the SOPT self-energy, eq.(^), and modify it into 
the form of eq.(^j) so as to interpolate between the weak and the strong correlation limits. Noting 
that Y^t(=t m )( n () = — we can easily determine the parameters A and B as 

A= n f [l-(N-l)n f ]^ 
»W - "/) 
1-2% 

U(N-l)n f (l-n f y ( 5) 

by the condition that S/(e) reproduces the correct high energy limit and the atomic limit, respec- 
tively. (We have used the relation Ef — Ef = U(N — l)n/.) 

The procedure proposed in this Letter can be summarised as follows. We start from the 
Hartree-Fock solution to G/(e) and calculate the second-order self-energy eq.(Q). The Fourier 
transformation method || is used in the calculation. We then construct the modified self-energy 
by eq.(||) with the Hartree-Fock term (N — VjUnj instead oiUrif. A and B are given by the above 
equations ( |l4| ) and (^), and rif — rif at this stage. To fulfil the Luttinger sum rule approximately, 
we subtract the constant term from the self-energy as eq. (|i"T|) . Next, we reconstruct G/(e) by 
eq.(|3|) using this modified self-energy, and calculate the cavity Green's function G/(e) by eq.(||). 
The self-energy eq.(f|) is calculated again by using G/(e), and repeat the above cicle until the 
convergence of the solution is achieved, rif and rif are recalculated at each iteration. 

As a numerical example, we show in Fig. 2 the f-electron density of states for N = 6, W = 1, 
V = 0.2, E = 0, U = 1 and (a) E f = -0.3 and (b) Ef = -0.5. The insets enlarge the low 
energy regions. The case (a) is metallic, while the case (b) is just on the border between the 
metallic and the (Kondo) insulating states. The number of f-electrons is (a) rif = 0.102 and (b) 
rif = 0.127, both of which satisfy the Luttinger sum rule within 3% of accuracy. Because of the 
strong correlation, the Kondo peak is renormalised to a small width and is further split by the 
mixing. One can see the large peak above the chemical potential, which represents the spectra for 
the transition f 1 — * f 2 and has a large weight due to the degeneracy. This peak and the peak below 
the chemical potential corresponding to the f 1 — * f° transition is absent in the SCSOPT calculation 
shown in Fig. 1 and can not be reproduced in the stronger correlation within that scheme. 

The new scheme proposed in this Letter can describe the electronic states of the strongly 
correlated electron systems with orbital degeneracy in a manner much simpler than that by Yeyati, 
et al.|| It can be easily combined with the density of states obtained from the band calculations, 
and will be used for an analysis of the realistic models of heavy Fermion materials with the strong 
correlation. 

The present scheme can be applied only to zero temperature when the Luttinger sum rule is 
forced to be fulfilled strictly. Since it is fulfilled only approximately in the present calculation, 
it is possible to apply the present scheme to finite temperatures to obtain semiquantitative but 
reasonable results. ||l2f 

Corrections of the next leading order of the lAi-expansion can be taken into account as a spin- 
fluctuation by using the technique proposed in |15j . Extension to the case with the crystalline field 
splitting, and calculations of the magnetic susceptibility and other physical quantities will be the 
important issues and are now under study. 
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Figure 1: The f-electron density of states for N — 6, W = 1 (the unit of energy), V = 0.25, 
Ef — 0.15 and U = 0.3 is shown. The straight and dashed lines indicate the SCSOPT and the 
Hartree-Fock results, respectively. 




Figure 2: The f-electron densities of states for N = 6, W = 1, V = 0.2, U = 
and (b) Ef = —0.5 are shown. The insets display the low energy parts. 
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